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ABSTRACT 



The canonical structure of pure Yang-Mills theory is analysed in the case when Gauss' law 
O ■ is satisfied identically by construction. It is shown that the theory has a canonical structure 
in this case, provided one uses a special gauge condition, which is a natural generalisation 
of the Coulomb gauge condition of electrodynamics. The emergence of a canonical structure 
depends critically also on the boundary conditions used for the relevant field variables. Possible 
boundary conditions are analysed in detail. A comparison of the present formulation in the 
generalised Coulomb gauge with the well known Weyl gauge (Aq = 0) formulation is made. It 
appears that the Hamiltonians in these two formulations differ from one another in a non-trivial 
way. It is still an open question whether these differences give rise to truly different structures 
upon quantisation. An extension of the formalism to include coupling to fermionic fields is 
briefly discussed. 
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1 Introduction 

This paper is concerned with analysing the canonical structure of pure Yang-Mills theory |l[ 
along the lines of two previous publications @, ||, in which a new gauge condition, called 
the generalized Coulomb gauge condition was introduced and used to obtain a straightforward 
canonical formulation of Yang-Mills theory, in the case when Gauss law is satisfied identically 
by construction. In the previous papers certain assumptions were made concerning the spa- 
tial boundary conditions of the Yang-Mills potentials. The boundary conditions in question 
are very important for the elucidation of the canonical structure. Here I will analyse the req- 
uisite boundary conditions in detail, and show that there is a set of boundary conditions at 
spatial infinity, which is consistent with the canonical structure. Rigorous proofs and heavy 
mathematical machinery are omitted here in the interest of simplicity. 

The basic quantity of Yang- Mills theory is the action S, which is given in terms of the gauge 
field G^y as follows, 

S = -\ Jd'x(G flu (A),G^(A)). (1) 

where the inner product ( , ) is defined below. For the gauge field G^ v I use a matrix notation 
(summation over repeated indices), 

G^{A) = GJ(A)T a = d u A^x) - VMx) - ig^x), A„(x)], (2) 

where the quantities T a are matrices in a convenient representation of the Lie algebra of the 
gauge group G, 

[T a ,T b ]=tf ab c T c , (3) 

and the quantities A^ are the gauge potential components, 

a,(x) = A;r a . (4) 

It is assumed that the gauge group G is semisimple and compact. The inner product ( , ) 
for any two Lie algebra valued (matrix valued) quantities A = A a T a and B = B a T a is then 
expressed with the aid of the Lie algebra structure constants f ab c as follows, 

(A,B) = h ab A a B b , (5) 

where 

Kb = -fab'' fbc' ■ ( 6 ) 

The quantity h ab and its inverse h ab are used to lower and rise Lie algebra indices, respectively. 
The notation used here is otherwise pretty standard or self-explanatory, with e.g. Greek letters 
used as indices denoting Minkowski space indices ranging from to 3, and latin indices from 
the middle of the alphabet (k, £,...) denoting space indices ranging from 1 to 3. The Minkowski 
space metric is taken to be diagonal, with signature (+,—,—,—). Unless otherwise stated, 
repeated indices are always summed over, be they Lie algebra-, spacetime- or space indices. 



It is convenient for future reference to write the action ([I]) in terms of a Lagrangian L. The 
action S is the integral of the Lagrangian L in an appropriate time interval [x®, x *]. Thus, 

S = |j ? dx°L, (7) 

where 

L = -\j v d 3 x (Gok(A), G ok (A)) -- A J v rf 3 x (G U (A), G m (A)) . (8) 

In the expression (Q) for the Lagrangian L, the quantity V is some appropriate domain in R 3 , 
which yet has to be specified. 

As is well known, requiring the action (|I|) to be stationary with respect to local variations 
of all the potential components A^, considered as independent quantities, yields the follwing 
field equations, 

V V (A)GT{A) = 0,^ = 0,1,2,3. (9) 

The "covariant gradient" V M (A) used above in Eq. (Q) is a convenient notion, 

V,(A) =8^ + 10^, ], (10) 

which will be frequently used in what follows. 

The non-Abelian Gauss law is obtained from the equations (|9|) for fj, = 0. Expressed in 
terms of the potential A the non-Abelian Gauss law reads as follows, 

V k (A)V k (A)A° - V k (A)A k = 0, (11) 

where 

A k (x) = d A k (x). (12) 

The time derivative of any quantity will frequently in what follows be denoted by a dot on 
top of that quantity, as in the equation (0) above. 

1.1 The Hamiltonian or Weyl gauge formalism 

It is well known J| that Yang-Mills theory can be expressed in a canonical form in the Hamil- 
tonian gauge, or so-called Weyl gauge | 4 = 0. 

The Yang-Mills Lagrangian in the case A = 0, which will be called Lyy here, is obtained 
from the expression (pi) by putting A = in that expression, 



i J d 3 *(A k , A k )-\J d 3 x(G M , G M ). (13) 



1 r „ , s :u. 1 
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The Lagrangian (|T3| ) describes a theory which is not quite the same as Yang-Mills theory, 
since Gauss' law is absent from the field equations following from the action principle with the 



expression ([13]) as Lagrangian. Gauss' law takes the following form in the case when Aq = 0, 
as seen from Eq. (|TTD above, 

V k (A)A k = 0. (14) 

The usual way to analyse Yang-Mills theory in the Weyl gauge A = 0, is to proceed from 
the Lagrangian (|T3|), disregarding Gauss' law to begin with. Using the variables A% and A\ as 
generalised coordinates and velocities, respectively, it is then perfectly simple to derive a canon- 
ical Hamiltonian formalism for the system defined by the Lagrangian (|i~3|) . The corresponding 
Hamiltonian H\y is, 

H w = -\J rf 3 x(7r fc , n k ) + l - J d 3 *(G ki , G ke ), (15) 

i.e. a simple sum of a kinetic term, depending on canonical momenta n k only, and a potential, 
or interaction term, depending on conjugate canonical coordinates A^ only. 



2 First attempt at a canonical formulation when Gauss' 
law is in force 

Gauss' law, Eq. (|TT|) above, will now be considered as an equation determining the (matrix 
valued) potential component Aq, for given space components A and A. This equation is a 
system of linear, elliptic partial differential equations with time x° acting as a parameter in an 
appropriate interval, with x G V being the independet variables. I will discuss the solvability 
of Eq. ([Ll]) subsequently, but proceed now by assuming the existence of unique solution A , 
which is a functional of the space components A and their time derivatives <9oA = A, i.e. 

A = A {a,a}. (16) 

The question is then whether the Yang-Mills system, which is originally defined by the 
action ([l]), permits a canonical structure when the potential component Aq is a solution to 
Gauss' law (|TTD, i.e. when Aq is given in terms of A and A by the expression (fCBD- It is 



possible to get some insight into this question without specifying the actual functional form of 
the relation ([16]) in minute detail. 

The Lagrangian of the Yang-Mills system, when the potential component A , is given by the 
relation (]16|) above, is obtained simply by inserting the solution ([16]) for A Q into the Lagrangian 
(|8]) above. The resulting Lagrangian will be called Lq, and is explicitly given as follows, 

L = ~J v d 3 x(v k (A)A {A,A}-A k ,V k (A)A°{A,A}-A k ) (17) 

-^d 3 x(G M (A),G M (A)). 

At this point it is appropriate to check whether the action principle involving the Lagrangian 
L in ([17]) above reproduces the field equations (P). It is perfectly straightforward to verify the 



following result, 

6 r dx°L = - p dx° f d 3 x (8A k , V (A)(V k (A)A° {A, A} - A fc ) - V e G M (A)) (18) 

i i 

- FJ dx ° L dV * ( SA ° { A - A ) • vk(A)A " { A - A ) - A *) ■ 

i 

Now the boundary conditions for the solution Aq to Gauss' law, i.e. the system of linear 
elliptic partial differential equations ([0]), enter into the discussion. If the surface term in 
Eq. fll8|) is non-vanishing, and not by itself a variation of some surface functional, then the 
Lagrangian ( |T7| ) is not a valid Lagrangian in the action principle which is supposed to lead to 
the field equations (§) for fj, = 1, 2, 3, when A is given by (|T6|) . 



In the first place the domain V is actually considered to be all of R 3 . The integrals over V 
will be given the following interpretation, 

d 3 x • •• = lim / d 3 x • • • (19) 

V ^^oo7| x |<R 



The vanishing of the surface term in Eq. (J18|) is then equivalent to the follwing, 

lim / dQR 2 (5A°[A, A], V {r) (A)A° f A, A} - A^) = 0, (20) 

R^°°J\x.\=R V L > ' 

where the superscript (r) denotes the radial component of the corresponding quantity. 



Thus, if the surface term (p0|) vanishes for all admissible variations of the independent 
generalised coordinates A and velocities A, respectively, then the variational principle 



d 



r x° 

/ f dx°L = 0, (21) 

Jx9 



leads to the following equations of motion, 

V {A){V k {A)A° {A, A} - A fe ) - ViG M {A) = 0, (22) 



as is evident from the relation (|T8|). Needless to say, the equations fl22|) are nothing but the 
field equations @ for fj, — 1, 2, 3, with A given by the formal solution fllED to Gauss' law (pT|). 



The vanishing of the surface term ( |20|) depends on the assumed asymptotic behaviour of the 
independent variables A and A, as well as on the boundary conditions (at spatial infinity) of 
the dependent variable Aq. I will return to this question below, and continue for the time being 



by assuming that the relation (2C) is valid. Then the Lagrangian Lq given in (17) ought to be 
a suitable starting point for the construction of a Hamiltonian and the corresponding canonical 
variables by means of a Legendre transform in the usual way. 

The formal definition of the canonical momentum P% conjugate to the coordinate A k a is, 

P fe V,x) = 6 J^ = (v fe (A)A {A,A}) V,x) -At(x°,x), (23) 



where the condition ( pH|) has been used in the calculation of the functional derivative of L$ 
in ( p3|) above. In view of the fact that A in Eq. (^) satisfies Gauss' law ([11]), one finds 



immediately from (|23|) that 

V fe (A)P fc (x°,x) = 0. (24) 

Now one is supposed to be able to solve Eq. ([£3) for the generalised velocity A% in terms of A 
and P, respectively. But this is impossible, since, Eq. (^) can not be solved for the quantity 
T defined below, 

r(x°,x) = V fc (A)i fc (x°,x), (25) 

i.e. if one tries to derive an equation for the quantity T defined above from Eq. (|23D , one gets 
a completely vacuous identity for this quantity, as a result of Eq. (|24|) . 



It would seem then, that the canonical formalism breaks down for the case at hand. However 
this is not necessarily the case. The difficulty described above can be avoided if one can manage 
to make Aq independent of the generalized velocity variables A. This can be accomplished in 
the present situation (Aq ^ 0) by imposing the following gauge condition, 

V k (A)A k (x°,x) = 0. (26) 

The condition fl2"6"|) is the generalized Coulomb gauge condition referred to previously. The fact 
that this is actually a proper gauge condition, has been demonstrated in Ref . . 

However, if one uses the gauge condition (|26|) then the generalised velocities A k are no 
longer independent quantities, and then one cannot use the formula (^3|) as it stands for the 
construction of the canonical momentum variables. An alternative procedure which leads to a 
proper canonical formalism will be given below. 



3 Gauss' law and asymptotic conditions 

I now assume that the generalized Coulomb gauge condition ( ^B] is in force. Then Gauss' law, 
Eq. (|ll|), takes the following form, 

V k (A)V k (A)A° = 0. (27) 

If one demands that the solution to Eq. fl2?D vanishes roughlty speaking faster than | x | _ 2 for 
| x |— > oo, then the solution vanishes identically. This can be seen as follows. Take the inner 
product of Eq. ( $H\ ) with A and integrate over x. Using the ordinary divergence theorem one 
readily obtains the following result, 

lim / d 3 x(V k (A)A ,V k (A)A ) = lim / dQR 2 (A , V {r) (A)A ). (28) 

Assuming now 

A)(*°,x) ||x[=* ~ R-" , V w (^)^o(^°,x) |, xM ~ R-^ 1 , 7 > i (29) 



one finds that the limiting value of the right hand side of Eq. fl28|) is zero. Since the inner 



product ( , ) is positive definite, one then concludes from Eq. (|28| ) that A is covariantly 
constant. However, since A vanishes at infinity by assumption, the covariant constant is zero, 
i.e. that 

v4 (x°,x) = 0. (30) 

Needless to say, the asymptotic condition (|29|) above can be somewhat refined; all that is needed 
to obtain the result ([$(]) is that the right hand side of Eq. (^8|) vanishes, which guarantees that 
Aq is covariantly constant, and then that Aq has the limiting value at infinity (or at some 
finite point), so that the covariant constant in question actually vanishes. 



However, the class of functions with (roughly speaking) the asymptotic behaviour (|2D| ) does 
not exhaust the class of possible solutions to Eq. ([27]). It is also possible to consider functions 
Aq which approach a non- vanishing constant matrix at space infinity, 

lim A (x°,x) = A = A a T a , (31) 

|x|— >oo 



where the real quantities A a are absolute constants. In addition to Eq. (plf) one can impose an 
asymptotic condition on the radial derivative of Aq, 

d 
lim | x | — — rAo{x°,x) = 0. (32) 

|x|^oo O X 



Requiring the validity of the asymptotic conditions ([H]) and fl32| ) and using standard arguments 
of potential theory || , one now derives an integral representation for the function A involving 
the (ordinary) Laplacian of that function, 

A)(x°,x) = A - -jW d 3 y . l y 2 y A (x°,y), (33) 

47T JR. 3 | x — y | J 

Using the the present form of Gauss' law, Eq. fl27D , one then converts ( |33"1) into an integral 
equation for the determination of A . For this purpose it is convenient to introduce some new 
notation, 

U ka b (x ,y):=2gf\ c A k %x°,y), (34) 

and 

V\(x ,y):=2gf\ c ^A kc (x } y)+g 2 r c , d r' ba ,A k d (x ,y)A ka '(x°,y). (35) 



d_ 
One then obtains the following integral equations, 



Al{x\ x) = A* - -L / r3 d'y^^yj \u k \{x\ y) ^|^ + V\(x°, y)A b (x°, y) 



(36) 



The integral equations fl36| ) constitute the starting point for the proof of existence of solutions to 
the present form of Gauss' law, Eq. (|27|). For this one needs naturally also to specify conditions 
on the potential components A%, which, together with their space derivatives determine the 
(unique) solution to Eq. (|56|). All this is a part of the "heavy mechinery" mentioned in the 
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Introduction, which I will omit in this paper. However it is appropriate to note the following 
asymptotic conditions, which are needed for the existence of solutions Aq to Eq. fl36|), 



A W> x ) - j^ > ^(*°< x ) - j^I > e > °- ( 37 ) 

I will now denote the solution of the system of integral equations (p^) by A {A}, and the 
corresponding Lagrangian by L 0Q (compare with Eq. (|i~7|) ) , 



J °° ~ 2 



l - J d 3 x (V fc (A)A {A} - A k , W k (A)A° {A} - A fe ) (31 



■\j v d'x((m 



d 3 x(G M (A),G kl (A' 



The Lagrangian fl3"8|) will now in the next section be used to derive the Hamiltonian for the 
Yang-Mills system under the condition that the generalized Coulomb gauge ( fffi) is in force. I 
implement this condition as a constraint, by means of a (matrix valued) Lagrange multiplier 
field C(x), which is used to modify the Lagrangian (j38| ) as follows, 

Loo -> £' = Loo + f rf 3 x(C(x), V k (A)A k ). (39) 

4 Canonical coordinates, momenta and Hamiltonian 

I now make a direct transition to a Hamiltonian formulation using the modified Lagrangian 
( |39| ) above, in a manner described in the general case by Berezin [0. The starting point is the 
familiar definition of canonical momentum variables n k , 

<(^ x) = I = (V k (A)A {A}) a -A a k - (V k (A)C) a , (40) 

0A a {X , Xj 

The equations ( ^Q|) above, together with the constraint equations (|26f) are now supposed to be 
solved for the quantities A k and C a in terms of the canonical coordinates A k a and momenta n k , 
respectively. Using Eqns. (26) and (|27|) , one finds immediately from Eq. ( |40D that 



-V k (A)V k (A)C = V k (A)ir k , (41) 

which, together with appropriate boundary conditions, defines the quantity C as an x-dependent 
functional of A and n, respectively, 

C = C{A,n}(x°,x). (42) 

It is certainly desirable that the solution C {A, n} to Eq. ([H]) be unique. The uniqueness is 
guaranteed if one uses the following boundary condition, 

lim / dQR 2 {C{x),V {r) {A)C{x)) = 0. (43) 

R^oo J|xl=i? 



The proof of uniqueness of the solution to Eq. ( f4ip under the condition (|^) has essentially 
already been given above in connection with Eq. (]27|). Namely, if there are two distinct solutions 
to Eq. (f4lD , then their difference satisfies the corresponding homogeneous equation, which is 
precisely of the form fl2T|). However, under the condition (^), the homogeneous equation in 
question has only the trivial zero solution, as demonstrated in the discussion following Eq. ( f27|) . 
Hence the solution C {A, -n} is unique. 

One now straightforwardly expresses the generalised velocity in terms of coordinate- and 
momentum variables, 

A% = (V fc (4) {A} - C)) a - <. (44) 

The construction of the Hamiltonian H then proceeds in the usual way. The relation defining 
the Hamiltonian H is the following, 

H= f d 3 x(7r fe ,i fe )-L o, (45) 

Jv 

where the quantity A ocurring in the expressions in (fS|) should be given in terms of canonical 
variables by the expression ((44]). It should be observed, that it is indeed the Lagrangian L 00 



which enters in the definition of the Hamiltonian H above, since at this stage the constraint 
is an identity. 



By straightforward calculation one finally obtains the Hamiltonian expressed in terms of 
canonical variables from the definition (]45|), 

H = -^J v d 3 x(ir k ,7r k ) + ^J v d 3 x(G M {A),G kl {A)) (46) 

d 3 x fa, V k (A)A° {A}) + I J v d 3 x (v k (A)C, V k (A)C) , 



!\ 



where the quantity C is the appropriate solution to the system of linear elliptic partial differ- 
ential equations (0), as discussed previously. 



5 The Hamiltonian equations of motion 

The Hamiltonian equations of motion are obtained by functional differentiation of the Hamil- 
tonian ( |46f) with respect to the canonical momenta and coordinates, respectively. 

The equations of motion for the coordinates are as follows, 

^ X °' X ) = TTT^a = (V*(A)(4, {A} - C)f - K, (47) 

which agree precisely with the expressions (f0|) as they should. In the calculation leading to 
Eq. ([17]) one encounters the following surface term, 

lim / d 3 xd k {C,V k {A)5 7r C) = lim f dQR 2 {C,V {r) {A)5 n C), (48) 

-R^°oJ|x|<K R-^ooJ 



which must vanish in order that the functional derivative in Eq. (|4"?|) be well defined. The 
surface term fl48|) in question vanishes as it should, in view of the general boundary conditions 
([43|) imposed on the quantity C {A, n}. 



The calculation of the functional derivative of the Hamiltonian fl4"6| ) with respect to the 
generalized coordinate variables A% is fairly lengthy, but nevertheless straightforward. In the 
course of the calculation one finds that a surface term analoguous to Eq. (JZQ) has to vanish in 
order that the functional derivative in question be well defined, i.e. that 

lim / d£lR 2 (5A°[A}M r) (A)A° {A} - A (r) ) = 0. (49) 

The condition ( f4"9"D is quite non-trivial. It should be remembered, that the variation of the 
quantity A Q can not be declared to vanish outside some finite region, since A is a dependent 
variable, determined by Gauss' law, i.e. in the present case by Eq. fl3"6|). Under essentially the 
conditions (|37|) with e > 0, I have been able to prove, that the iterative solution of the Eqns. 



(36), which are equivalent to the system of partial differential equations (Q7u with the boundary 



conditions fl3T|) and fl32|), is such that, for large | x 

M (A)(x°,x) = O(^], (50) 




for any local variations 5 A. If one assumes that e > 1 in the discussion above, it is essentially 
trivial to show conclusively that Eq. ( |50f ) is valid. However, I will proceed by assuming that 
Eq. ( |50"D is valid also if one merely takes e > 0, which is a plausible assumption as discussed 
above. 



If the conditions fl3~7|) and (|50|) are valid, or more precisely if the condition ( (49"|) is in force, 
then one obtains, 

r tt 

n k (x°,x) = - - = - l g[A {A}-C 1 V k (A)C + 7r k } a (51) 

SA%(a ; u ,x) 



V e G H {A)) -tg[C,V k (A)A {A}] a . 

/a 



The pairs of equations, (|47D and (|5T|) are supposed to be equivalent to the original field 
equations ([5J). I will analyse this equivalence below. For this purpose it is convenient to note 
that the equations fl4"7|) and (|5T|) admit a constant of motion, 

d (V k (A)n k ) = 0, (52) 

which is crucial for the establishing of the equivalence between Eqns. (|9]) and the Hamiltonian 
equations of motion (f47D and (|5lD . Recalling the equation ([|l]) defining the quantity C together 
with the boundary conditions (0), one finds using Eq. (|52| ) that 



V k (A)V k (A)C = K 1 , (53) 



where K x denotes a constant. 



So, if the equations of motion for the canonical variables A and n, respectively, are in 
force, then the quantity C is determined by the equation (0), where the constant K\ is so far 
unknown. 

Eliminating the canonical momentum variables ti between the equations (HTj) and (EH]) one 
obtains the following equations, 



\k 



V u (A)G ku (A) + V k (A)C + 2ig \A\ C\ = , k = 1, 2, 3. (54) 

At this point one should recall that Gauss' law is in force by construction, 

V k (A)G ok (A) = 0. (55) 

Were it not for the terms involving the quantity C in Eqns. ([54]) one could now conclude that 
the Hamiltonian equations of motion ( f4"T| ) and ( pT| ) are equivalent to the original field equations 



(H). In order to obtain the required equivalence one must demand that 

V k (A)C + 2ig[A k ,C] =0, (56) 

simultaneously with the Hamiltonian equations of motion. In this situation the quantity C is 
determined by the partial differential equations (p3|). However, it does not appear to be possible 
to prove that the unique solution to (|53|) also satisfies (|56|) if the constant Ki has an arbitrary 
value, different from zero. However if 

Ki = 0, (57) 



then substituting this value in Eq. ( p3|) one obtains the final equation determining the quantity 
C, when all the other equations of motion are in force, 

V k (A)V k (A)C = 0. (58) 

But Eq. ([)8|) has only the trivial solution C = in the class of functions satisfying the boundary 
conditions ( f£3| ) as discussed previously. Hence if the constant of integration K\ in Eq. ( |5*B|) 
equals zero, then Eq. ( [56"D is trivially true, whence the Hamiltonian equations (f47D and (|51~[) are 
completely equivalent to the original field equations (Q). Thus by demanding that the condition 
( |57|) be valid, one obtains complete equivalence between the original field equations @ and the 
Hamiltonian equations of motion ([|7]) and (|5lD which have been obtained using the generalized 
Coulomb gauge condition (|26|). 



It is now appropriate to return to the Hamiltonian ([46]). This Hamiltonian differs in a non- 
trivial way from the Hamiltonian in the Weyl gauge, i.e. from the expression (|i~5"D , mainly due 
to the C -dependent terms in the former expression, but also due to the fact that a non-trivial 
A -dependence is possible in the case of the Hamiltonian (flop . The fact that the functional 
C {A, it} actually becomes zero when all the equations of motion are in force, does not mean 



that one can set C equal to zero in the Hamiltonian (|46|) and use the resulting expression 



to generate the appropriate equations of motion by means of functional differentiation. The 



functional dependence of the quantity C { A, ir} in the Hamiltonian ( f46"D is essential in order to 
generate the proper equations of motion. 

10 



Finally I comment briefly on the possibilities to generalize the considerations in this paper 
to a situation in which one couples the Yang-Mills field to e.g. a fermionic field. The field 
equations @ then get replaced by the following, 

V V {A)G>»{A) = r = g^T a ^T\ (59) 

to which one has to add the field equations for the fermionic fields. Also in this case is it 
possible to use the generalized Coulomb gauge condition ( |26| ) so that the Gauss' law takes the 
form 

V k (A)\7 k (A)A = J . (60) 

The construction of the canonical variables and Hamiltonian in this case, along similar lines 
to those presented here for the pure Yang-Mills case, does not meet with any difficulties of 
principle. Likewise, coupling the Yang-Mills field to a scalar field also gives rise to a system 
which has a canonical structure. Details of these constructs as well as some mathematical detail 
only briefly touched upon in this paper will be given in future publications ||. 
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